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Abstract 



The problem of splitting effects by vertex angles is discussed for non- 
integrable rational polygonal billiards. A statistical analysis of the decay 
dynamics in weakly open polygons is given through the orbit survival prob- 
ability. Two distinct channels for the late-time relaxation of type t^ are 
established. The primary channel, associated with the universal relaxation 
of "regular" orbits, with 5 = 1, is common for both the closed and open, 
chaotic and nonchaotic billiards. The secondary relaxation channel, with 
(5 > 1, is originated from "irregular" orbits and is due to the rationality of 
vertices. 
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I. INTRODUCTION 



Polygonal classical billiards is an active subject of research in mathematics and physics 
(see for review Ref. [|l|). In view of the null Lyapunov exponent and the null Kolmogorov 
metric entropy the rational polygons, formed by the piece-line billiard boundary with the 
vertex angles are rational multiplies of vr, are known to be nonchaotic systems |T]-§[. They 
are therefore well distinct from the Sinai billiard (SB) and the Bunimovich billiard Q 
(BB) where classical chaotic motion regimes are due to, respectively, dispersive effects caused 
by the circle disk and the squared boundary, and the interplay between boundary segments 
formed by the circle and the square. Meanwhile, the rational polygons of m equal sides 
and equal vertices (hereafter, the m-gons ^) revealed positive Lyapunov exponents with 
increasing of m. Furthermore, the polygonal billiards exposed chaoticlike changes in the 
associated quantum-level spectra [T^, which fluctuations are shown [|ll|)|ll] to be very close 
to the Gaussian orthogonal-ensemble-type statistics. In view of the splitting effects by the 
angle vertices, polygons do not satisfy the conditions of integrability PJT3|1 

These evidences for the chaoticlike features of the nonintegrable rational polygons were 
recently questioned by Mantica through the orbital complexity analysis. Unlike the case 
of the integrable billiards [0 , a delicate problem of the interplay between the regular (piece- 
line) and the irregular (vertex-angle) boundary segments in polygons cannot be solved in 
terms of the first-order averaged polygonal characteristics, such as the average-orbit cod- 
ing length 1^, the mean wall-collision time [|14|, or the average collision number [|l^]. As 
shown through analysis of the orbit- wall-collision statistics [|^, the higher order correlation 
effects induced by vertices play a crucial role in the intrinsic dynamics of m-gons. Earlier 



this was experimentally corroborated |Tl[ through the quantum-level spectra statistics: "un- 
like the case of the irrational polygons, the long-range level correlation effects are due to 
the rationality of vertex angles" In the case of the curved-by-circle infinite-m ratio- 



nal polygon (oo-gon) the late-time memory effects arise from the sliding orbits [|T5[, which 
have no analog in the ballistic-type dynamics of its counterpart given by the circle billiard 
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(hereafter CB), where the orbit classification is well established [|^. This implies that the 
quasi-classical approach has no justification for the "quantized" vertex-spliting oo-gon that 
is geometrically equivalent to the CB. In other words, the quantum-to-classical dynamics 
correspondence suggested between a given m-gon and its circumscribing counterpart is vi- 
olated |T5| regardless of the fact that the geometric correspondence exists and can be achieved 
with any precision when m — oo (with the help of the aforementioned first-order averaged 
characteristics) . This is in agreement with a conclusion on inapplicability in polygons of the 
quantum-to-classic correspondence principle. The latter was elaborated within the scope 
of the conventional Wentsell-Kramer-Brillouin picture that failed to establish a one-to-one 
correspondence between classical orbits and their quantum counterparts. 

The issue of the current paper is an investigation of the vertex-splitting effects revealed 
in the decay dynamics of finite-m rational polygons. An analysis is given through the sur- 
vival orbit statistics in the weakly open billiards, which boundaries permit orbits to escape 
through a small opening. We will see that, similarly to the case of the intrinsic dynam- 
ics in the closed m-gons |T3, the vertex-splitting effects are dual of vertex-ordering and 
vertex-disordering motion effects, that is manifested through, respectively, the orderlike and 
chaoticlike behavior observed in the late-time relaxation. The paper is organized as follows. 
The recent findings for the decay dynamics in nonchaotic [l^ and chaotic [|l^] billiards are 
given in Sec. II and discussed within concepts of the primary and the secondary relaxation 
channels. The weakly open rational polygons are analyzed numerically and analytically in 
Sec. Ill for the cases of small and large numbers m. Discussion and conclusions are summa- 
rized in Sec. IV. 



II. DECAY DYNAMICS OF THE CHAOTIC AND NONCHAOTIC BILLIARDS 

The intrinsic dynamics of the closed classical billiards is commonly considered in terms 
of a temporal decay of the correlation functions for certain dynamical variables (see e.g. 
Ref. []TB|). Pure exponential lost of amount of memory on the initial states is not a unique 
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channel of relaxation even in the chaotic systems (see e.g. Ref. [|T9[)- By studying the chaotic 



billiards, such as the SB [|T^-p2[ , that is dynamically equivalent to the correspondent Lorentz 



gas (LG) model |]I8[ and the BB P J23| , |2^ with a stadium geometry, it has been recognized 
that a crossover from the short-time exponential to the late-time algebraic decay is due to 
the long-term memory effects on a free regular motion. The algebraic tail of the correlation 
functions seems to be vanished only in the case of the fully hyperbolic systems correspondent 
to such geometries as the finite-horizon SB 0,^] (equivalent to the high-density LG) or 



the diamond . Qualitatively the same can be referred to the decay dynamics in weakly 
open billiards that describes a crossover from a bounded to unbounded orbit free motion. 
Such a decay dynamics is initially established by uniformly distributed Nq point particles 
(of unit mass and unit velocity) moving inside of the closed planar billiard table and then 
allowing to escape through a small opening of width A. A temporal behavior of the dynamic 
observables can be scaled to the characteristic billiard times, namely 

T"c = -p- and Te = —Tc. (1) 



Here the mean collision time [ P^p5| , pB[ Tc and by the mean escape time |T^,|TB|J2^,|2^ are 
introduced through the accessible area A and the perimeter P (^ A) for a given billiard 
table. The late-time (t ^ Tg) algebraic-type evolution of non-escaped orbits (or particles) 
Nt cc NQ^Te/ty, which time derivation is related to the survival probability, is definitely 
characterized by the decay dynamic exponent 6. 

In the integrable nonchaotic billiards, the first observation of the algebraic decay with 



exponents 5 ~ 1 was given for the case of the square billiard. Recent study of the decay 



dynamics in the 4-gon revealed [|Tj] two distinct channels of the algebraic-type slow relax- 
ation. The first one is due to the regular-orbit motion with the decay dynamic exponent 
6 = 1, and the second channel is originated from the "irregular" orbits induced by the sin- 



gular vertex-spliting effects, which give rise to the subdiffusion motion regime indicated ||n 
by 6 ~ 0.85. Such a kind of diverse dynamical behavior was observed through the survival 
orbit spectra defined by a number for the survived orbits and simulated in the integrable CB 



and in the almost-integrable square billiard (see, respectively, Figs. 5,6 and 3,4 in Ref. 
|T^). In both the cases an irregularlike motion is due to orbit families known as, respec- 
tively, the "whispering-gallery" (polygonlike) and the "bouncing-ball" orbits. Meanwhile, 
in the late-time decay of the CB (see inserts in Figs. 3,4 and 5,6 in Ref. [0) the short-time 
living, frequently escaped "whispering-gallery" orbits, unlike the nonintegrable case, did not 
contribute to the second channel of relaxation. Thus, the solely exponent 6=1 was observed 
in the CB. 

In the chaotic closed and weakly open (including Hamiltonian) classical systems, pre- 



sented by the BB p|,|4|,|8|,|9[ , the infinite-horizon SB |T3Jl|,|20|,||,|30|,|l , and by the low- 
by the 



density LG model pl] , |32| , |33| , the algebraic- type decay was numerically revealed 
dynamic exponents 6 > 1. Similarly to the nonchaotic case, it has been repeatedly rec- 
ognized that the algebraic tail is caused by the "arbitrary long segments" observed in the 



evolution of stochastic orbits ||2J], or by the regularlike orbit motion due to "sticking parti- 
cles" [^,0. This implies that in both the cases this relaxation is due to a free motion of 
the corresponding trajectories in the infinite (with respect to the observation time) distinct 
corridors which are open in the relevant phase space |]SD|-^ . 

The algebraic-type relaxation channel with 6 = 1 {= a) established in the chaotic 
||T7| , [l9| , p3| , |30|j3l| as well as in nonchaotic ll^Jl^ billiards, seems to be generic for all non-fuUy 
hyperbolic systems with smooth convex boundaries. Its independence of the billiard space 
dimension [|31[, its insensitiveness to details of the boundary shape [1^, including that to 
a position of the small opening p3|, and to the initial conditions [|14|, suggests that the 
late-time a-relaxation arise in classical systems as the universal primary relaxation. The 
latter is a part of the two-step relaxation scenario discussed in the chaotic [|17| and non- 



chaotic weekly open classical systems. This universal scenario was introduced by the 
short-time pure exponential and by the first-power-algebraic overall orbit decays given by 
Nt cc Ai'oe"*^*/^'^^^ with 7 = 1, and by Nt cc NQ{Te/tY with 5 = 1, respectively. Meanwhile 
the real relaxation does not exclude the existence of the additional intermediate transient 
regimes approximated by the stretched-exponential form with 7 < 1 (earlier discussed 



for chaotic billiards in analytical and numerical |T8| forms), and by other algebraic-decay 



forms pO| , |3T| with 5^1. The escape mechanism of the primary relaxation in the chaotic 
and nonchaotic billiards was described in details within a coarse-grained approximation (see, 



respectively, Eq.(12) in Ref. [|r^ and Eq.(15) in Ref. ||14[)- Unlike the case of the primary 
relaxation, the temporal observation conditions (observation windows) for the secondary re- 
laxation in the chaotic billiards {S = (3 > 1) are shown to be very sensitive to the billiard 
geometry [l^, to the dimension [^] d of a billiard table, and to the initial conditions [jl4 



This can be exemplified by the dynamic-exponent constraint 1 < f] < d proposed in Ref. 



| 31| and observed in the chaotic BB and SB |19| , |30| , |3lll billiards 



The algebraic-type decay of correlations in the low-density chaotic LG is due to evolution 



of trajectories within the infinite principal and/or "hindered" open corridors pO|-^. In the 
corresponding disk-dispersing SB (of side L and disk radius R) with the infinite-horizon 
geometry {R < L/2) the principal and the "hindered" corridors, respectively, can be gov- 
erned by the limiting disk radii Ra = L/2 and Rfs = ^/2L/4: and observed through the 



established |T^ algebraic-tail temporal windows. Thus, the primary relaxation was observed 
within the domain R13 < R < Ra, when the "hindered" corridors are closed Under 
the geometric conditions < R < Rj3 the transient /^-relaxation channel was activated and 



indicated by the decay dynamic exponents P = 1.2 and /? = 1.1, respectively, in Refs. ||30 
and 



One can see that the primary and the secondary relaxation channels can be geo- 
metrically associated with Bleher's principal and "hindered" corridors, respectively. The 
latter case was additionally characterized as a superdiffusive motion regime, observed 
earlier in the chaotic [p!8| , p0| , p2| , p3| and very recently in the polygonal [0 closed billiards. 

We see that non-fully hyperbolic billiards, on one hand, are indiscernible within decay 
dynamics observed through the primary universal a-relaxation channel. On the other hand, 
the chaotic and nonchaotic, the open and closed billiards are well distinguished with respect 
to the secondary /5-relaxation given by the decay dynamic exponents (3 > 1 and (3 < 1, 
respectively. In what follows we give numerical and analytical analysis of the conditions 
of stabilization of both the primary and secondary relaxation channels in the weakly open 
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rational polygonal billiards. 



III. ORBIT DECAY IN POLYGONS 



We deal with the rational polygons of m equal sides, denominated as m-gons, circum- 
scribed below a circle of radius R. The mean collision time r^m = (vr/2/2) cos(7r/m) and 
the mean escape time = (7ri?^m/2A) sin(27r/m) are given with the help of Eq.(|I]) 
through area Am = (mi?^/2) sin(2'7r/m) and perimeter Pm = 2mi? sin(7r/m). In the limit 
m (yo one naturally arrives at the circle geometry of the oo-gon with the mean times 
Tcoo = TcR^^ = 7ri?/2 and t^oo = T^'ji^^ = 7r^i?^/A , both characteristic of the CB. This 
demonstrates how the first-order dynamic characteristics can be introduced through the 
aforementioned geometrical correspondence that takes place between the oo-gon and the 
CB. Meanwhile the dynamical correspondence does not exist |jT5| because the vertex- memory 
effects violate commutation between of the temporal (t oo) and spatial (m oo) limits. 



A. Small Number of Vertices 



Similarly to the closed m-gons [T^, let us consider the case of small number of vertices, 



with m < 10 , within the scope of the deterministic approach. This is given by generalization 



of the regular-orbit description introduced [|14| for the particular case of m = 4 and is 



straightforwardly based on the fact that the wall-collision angles ip (counted off the normal 
to the boundary and preserved by elastic refections) are integrals of motion. This is true 
for the integrable billiards where dispersing or splitting effects are absent. A dynamic 
description of the regular-orbit motion can be introduced with accounting of the fact that 
m (or m/2) sides of a given m-gon, with odd (or even) number of vertices, are dynamically 
equivalent [0. A description of the wall-collision statistics can be therefore reduced to the 
collision- angle domain Lp = [0, ipm] with 
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7r/2m for odd m, 

= { (2) 

7r/m for even m. 

In turn, the yj-family regular-orbit sets can be introduced (for details see Appendix) through 
the characteristic collision times, namely 

TxR cosn/m simpm 

tcmm = 7; 7 , with (3) 

{0 for odd m and m/2, 
(4) 
vr/m for even m/2. 

The collision time tcm(v^) is related to the billiard mean collision time through the 
"mean-collision-time" equation, i.e., < tcm{f) >c= Jq"" tcm{f)fom{'^)dip = Tcm- The latter 
is equivalent to the known "mean-free-path" equation considered in the uniformly pop- 
ulated two-dimensional {2D) collision subspace Qcm (see also Eqs.(3,4) in Ref. the 3D 
phase space Q^- The y^-set-orbit distribution function /om is defined here as 



cos(y-t/'^) 



for subspace flc 



fornix) ={ , ; ^ (5) 



by generalization of Eqs.(6,7) in Ref. [|14 . 



We discuss the late-time [t ^ Tem) survival dynamics in a given m-gon through the cp-set- 
orbit decay spectra defined by numbers of the survived orbits Nfn{t, f) and by corresponding 
overall numbers Ntm =< Nm(t, if) >c- Here a procedure of averaging over collisions, denoted 
by < ... >c, is introduced through the aforesaid "mean-collision-time" equation. The univer- 
sal relaxation channel, associated solely with the regular orbits, is given through numbers 
of the nonescaped orbits, namely 

— = Cm[^) fomi'p)— and -— = D^—, (6) 

-''Om Tcm -''Om 

which are the late-time solutions of the relevant decay-kinetics master equation (see Eqs. 
(16,17) in Ref. JI^)- As seen from Eq.(^, the fundamental characteristics tcm(v^), fom{f) 
and are common for both the decay and intrinsic dynamics. The (p-set orbit-partial 
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weight Cm{,^) and the orbit-overall weight Dm of the corresponding algebraic tails can be 
established in explicit form within a certain coarse-grained scheme and directly observed 



l^ , p!7[| . However an interesting analysis that indicates a departure of m-gons from the true 
integrable systems due to vertices can be given through the algebraic-tail weights without 
detail calculations. 

By taking into account the mentioned equation Ntm =< Nm(t, ip) >c, the overall regular- 
orbit-decay weight Dcm = Ccnf' =< Cmi^p) >c^^'^ immediately follows from Eq.(||). The 
validity of this relation can be straightforwardly tested for the integrable CB where the 
irregular orbits do not survive in the late-time relaxation. Indeed, in this case C^^) = 0.206 
and -D^i?^^ ~ 0.214 that follows from the observation (see Tab. 2 in Ref. JI^) of the primary 
relaxation of the collision space QcR ■ One can see that the algebraic-tail parametrization 
of the regular-orbit decay effects in the CB is self-consistent that is experimentally justified 
with a high precision, i.e., C^^^^^ = -D^^^"* = 0.210 ±2%. In the case of polygons, a violation 
of the relations Ccm^"* = Ccm^^ and Ccm^^ = D^cm'^ is expected in view of the long-living 
irregular-orbit motion induced by vertex-splitting effects. 



We have performed numerical experiments on decay dynamics in m-gons with small 
number of vertices: m = 3, 4. ..8. The initial particles {Nq^ = 10^) have been distributed 
randomly within the two distinct phase spaces described in Eq. (^) and then allowed to escape 
through a small opening A ^ i?. The condition Tem (= 300) is chosen to be common for 
all m and that has been provided with the accuracy of ±5% with the help of Eq.([l|). In all 
cases of m < 8 the typical algebraic decay is observed within a certain temporal windows 
given by, approximately, lO^Tem < t < lO^rgm- The particular cases of the observed decay 
spectra for pentagon and heptagon are exemplified in Fig.l. In general, the overall-orbit 
late-time decay in m-gons with small number of vertices do not show noticeable deviations 
from the linear relaxation ||39| (see the left insert in Fig.l). Thus the partial weights Cm^^\^) 
are derived from the observed numbers Nm^^\t,ip) through Eq.(^ with accounting of the 
estimated distribution function fomi'^) and the collision times tcmif) given in, respectively, 
Eqs.dl) and (|). These equations were additionally experimentally tested {e.g. see the 
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right insert in Fig.l). 

As seen from Fig.l, the observed partial weights Cm^^\^) exhibit regular (small) and 
irregular (large) deviations from the mean magnitude Cim^^ = Cjm'^ shown by a solid 
(horizontal) line. The latter and the regular-orbit weights Cim^^ (evaluated with regardless 
of the large isolated peaks) are accumulated in Table. 1. Through analysis of a difference 
between the overall and partial weights ACcm = Cctu^ — Ccnf' ■, one can see with the help 
of Table 1 that, similarly to the case of the intrinsic dynamics (see Fig.2 in Ref. |]l5l), the 
vertex-splitting effects in the even-gons are more pronounced than those in the odd-gons. In 
all the m-gon cases a deviation of the total weights, treated as ADcm = Dcm^ — Dcrif* (with 
-Dcm*'* = Dim^^ and Dcm^^ = Ccm^^ given in Table 1), exceeds experimental error (±2%) 
established above for the integrable SB, and we therefore infer that Z)cm*^ > Dt^\ This 
implies that the irregular-orbit motion is involved into the observed relaxation. On the 
other hand, no noticeable deviation from the primary relaxation is indicated in the observed 
decay-orbit dynamics when m < 8. Similarly to the case of the intrinsic djTiamics observed 
PJT3[| in the m-gons with small number of vertices, we deduce that the regular-orbit motion 
dominates in the late-time decay dynamics. 



B. Large Number of Vertices 

The universal two-step relaxation in open billiards is shared by m-gons with arbitrary 
number of vertices |^0|- Figs. 2 and 3 we analyze our numerical results for the late-time 



overall-orbit decay dynamics in the m-gons with large number (m = 2" with n = 3, ...6) 
of vertices for the cases of relatively small and large opening widths. In general, one can 
observe that the decay dynamics of the rational polygons with increasing of number of sides 
moves away from that given by the geometrically corresponding CB: the chaotic effects, 
manifested by the secondary relaxation channel with (3 > 1, become more pronounced with 
number of vertices. As seen from the observed relaxation of the initially equivalent states 
(given in phase spaces floo and Qr) is quahtatively distinct, and the open oo-gon and the CB 
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are not therefore dynamically equivalent. Moreover, the m-gons with m > 8 do not expose 
the algebraic decay with /3 < 1, characteristic for subdijfusive motion regimes in nonchaotic 
open systems observed |]T^ in 4-gon. 

In the particular case of A = 0.05i? shown in Fig. , the universal relaxation remains 
stable until m = 64, but when m > ma^"^^ = 128 the primary relaxation channel turns 
up to be closed. In other words, the regular-orbit-motion relaxation affected by vertices is 
transformed into the irregularlike- motion relaxation indicated by dynamic exponent f3 = 1.2. 
Qualitatively the same follows from Fig. 3, but the upper limit for the a-channel observation 
window shows its dependence on A, i.e., rria^^^ = 32 for A = 0.20i?. Thus Bleher's principal 

corridors of the regular-motion relaxations are open for 3 < m < mo^^'' (similarly to the case 
of the SB given hy Rp < R < Ra that is discussed in Sec. II) when the vertex-splitting (or 
disk-dispersing effects) are weak |15]. We see that the a-relaxation ( 3 < m < m^^^'') occurs 
as a precursor of the /3-relaxation regime, realized for m > 'mt^^\ just as in the chaotic SB 



case the /^-relaxation (0 < i? < Rj^) was observed |40| before the a-relaxation [Rp < R < 
Ra)- Conversely, the observed a-to-/5-relaxation crossover in the m-gons, induced by the 
interplay between the piece-line regular and the vertex-angle singular boundaries, is similar 
to that between the semi-square and semi-circle parts of the stadium boundary. The latter 



was deduced in the case of the BB with small (A = 0.01) and large (A = 0.25) opening 
widths, with the observed exponents 1 ^ /? ^ 2 (see Table 1 in Ref. [^). 

Qualitatively, the effect of closing of the principal corridors in a given m-gon can be un- 
derstood by difficulties, increasing with m, to draw the long segments of free motion, which 
intersect the polygonal sides in the correspondent LG lattice but avoid the vertex angles (see 
also discussion in Appendix). By contrast, the /3-relaxation revealed in Figs. 2, 3 is asso- 
ciated with stabilization of the "irregular" -type trajectories, which are effectively deviated 
by vertex angles. More precisely, the observed order-to-chaoticlike crossover can be under- 
stood as a regular-to-irregular orbit transformation of the aforesaid sliding orbits (formed 
||r5| hj if ^ 7t/2 -sets with characteristic free-path times t^cm^\ip) = Tcoo cos~^ following 
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from Eq.(^), which survive in the open m-gons with m < rric into those, renormahzed sub- 
stantially by vertex angles, called |jT5[ by vortexlike orbits (with the large, but finite mean 
characteristic time Tcm^^^^ = Tcoom/n), which expected to be stable for m > nia ^ 1- Within 
this context, the survival conditions of the regular (sliding) and the irregular (vortexlike) 
orbits driven, respectively, by piece-line and vertex-angle parts of the open polygonal bound- 
ary, can be introduced as follows. On one hand, the favorable observation conditions for the 
a-relaxation (or the /3-relaxation) should ensure to exclude (or to include) the vertex-angle 
effects under the constraint m < rria (or m > rria). On the other hand, in the weakly open 
(A ^ Pm) m-gon of a side length Lm = Pm/^, the survival conditions for the regular sliding 
or irregular vortexlike orbits are satisfied by geometric constraints, respectively, A ^ 
or A ^ Lm- Hence, the a-/3-relaxation crossover, observed at m = rria, is ensured by the 
condition A = L^- With taking into account that perimeter in the polygons with large 
number of sides is well approximated by Pm ~ 2ttR, one arrives at the desirable criterium 

rria = (7) 

This finding provides the estimates rria^'^"'^'' = 126 and 31 for the particular cases of ma^"^^ = 
128 and 32 realized in Figs. 2 and 3, respectively. We infer that unlike the case of the closed 



polygons |15[, the vertex-splitting effects in the open m-gons with large number of vertices 
give rise to stabilization of the vortexlike-orbit motion. 



IV. DISCUSSION AND CONCLUSIONS 



The mild discontinues caused by vertex angles and relative lengths of the edges is the 
central problem of the intrinsic dynamics of the "almost integrable" polygonal billiards 
commonly discussed |^ in terms of the orbit ergodicity, mixing, entropy, coding, complexity 
0], pseudo-integrability 0, spectral level statistics pO|JTl| , and the orbit collision statistics 
p!5| . The problem is now addressed to the decay dynamics in the m-gons and is discussed 
through the orbit survival probabihty \l'm(^) = \d{Ntm/NQm) / dt\ related to the number of the 
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survived orbits Ntm- The decay spectra given by the y^-set regular-orbit numbers Nm(t,ip), 
are also studied for the case of small number of vertices m. 



A general approach to the decay problem based on a simple decay kinetic equation |T7 



naturally arrives at the primary slow relaxation of the regular-orbit sets given by a 
(see Eq.(^). We have demonstrated that the universal a-channel, attributed for both the 
chaotic [jl^ and nonchaotic billiards, is also characteristic of nonintegrable rational 



polygons. The primary relaxation-motion regime originated from the piece-line parts of 
the polygonal table is associated with long-living sliding orbits with large collision angles 
(p. In the corresponding phase space these orbits are unbounded trajectories (see Fig. 4) 
that move without splitting at angle vertices along Bleher's principal corridors. Following 
to the simplest polygonal orbit classification by Gutkin the regular sliding-orbit sets 
can be presented by the "infinite-past-to- infinite-future" trajectories. They "never" hit 
vertices, preserve the initial linear momenta, and show a regular behavior in the orbit decay 
spectra N^it, ip) (see Fig.l). Conversely, the singular orbit sets caused by the " infinite-past- 
to- vertex", the "vertex-to-infinite-future", and the diagonal "vertex-to- vertex" trajectories 
^ exhibit pronounced weights Cm{f) in the orbit-decay process (shown by high peaks in 



Fig.l). Eventually, they do not play any significant role in the wall-collision statistics in 
m-gons with small number of vertices limited by, approximately, 3 < m < 8, and thus the 
primary relaxation dominates. This corroborates by our numerical study (analyzed in Table 
land Fig.l) and, in general, is in accord with studies of the closed polygons by the orbit-wall 



collision statistics |I5| and by the orbit complexity [g]. 

When the number of vertices is large, the secondary relaxation with the survival probabil- 
ity \l/m^ oc t~^~^ predominates over the primary relaxation (see Figs. 2,3). The established 
domain for the decay exponent 1 < j3 < 2 corresponds to that known |^ for the chaotic SB. 
Qualitatively, the survival probability function can be associated with the distribution 
function for trajectories trapped by the strange attractors, discussed in the theory of the 
open classical chaotic systems, or with the corresponding waiting-time probability function 
43| . With accounting of findings for the decay dynamics on the SB by Fendrik's group 
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30| , |3H] , one can expect that the secondary relaxation is due to the singular trapped orbit 



sets that move freely along Bleher's "hidden corridors". Similarly to the case of the chaotic 
billiards, the observation conditions for the secondary relaxation in rational polygons are 
sensitive to the initial conditions and to the geometrical constrains. Indeed, the /3-relaxation 
channel turns up to be closed if the initial particle distribution is simulated |^ in the 2D 



collision subspace Qcm- In the case of the 3D Qm space the secondary relaxation appears 
to be dynamically stable under the geometrical constraint m > nia, where rria = 2'n'R/A 
is given by the a-to-/3-relaxation criterium estimated in Eq.(0). As shown, this criterium 
meets the favorable survival conditions for the regular-motion regime with those induced 
by rationality of the vortices. The latter are generated by the sliding orbits through the 
vertex-" ordering" effects and are associated with the vortexlike orbits [1^. As follows from 
Eq.(0), the observation window of such a motion disappears in the limit A — 0, when the 
vortexlike orbits do not survive in the closed polygons (see Fig. 2 in Ref. |jl5|]). Finally, we 
have demonstrated that the vortexlike orbits become stable in the open rational polygons 
and visible through the secondary slow relaxation common for the chaotic systems. 
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V. APPENDIX. ORBIT-SET COLLISION TIMES 

In a given m-gon a number of geometrically equivalent walls k is bounded above by 

{m for odd m; 
(Al) 
m/2 for even m. 

The current collision angle cpkm with a wall k (= 1, 2...qm) of a y^-orbit with = [0, 7r/2gm] is 
reduced through the relation ipkm = f — ^km with the help of 6^^ = [— 7r/2, 7r/2] introduced 
as the lowest angle between the /c-wall and axis x , namely 
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gm — 2A; + 1 for odd g^; 



0^™ = ^ < (A2) 

Q'm ^ 2 A; for even g.^. 

As shown in Fig.4 for the particular case m = 3, the estimates for the wall-colhsion times 
tcmif) are found through summation of numbers of intersections n{t,ipjsm) for a trajectory, 
induced by a given <^-set orbit, considered in the correspondent infinite LG lattice, namely 

qm 

The estimation procedure can be exemplified by a relation tcos(y9i3) = n{t,(fi3)3as. The 
latter employes the fact that a distance between the equivalent walls is 803, where = 
Rcos{7i/m) stands for the apothem in a given 3-gon. This yields 



ncm{t, ^) = rj—j;^ = "-(^5 ^km) (A3) 

k=l 



^ C0S(V3 - Qkm) 
.k=l 



(A4) 



where qm and 9^^ are given in Eqs.(Al) and (A2), respectively. Straightforward estimation 
of Eq.(A4) [0 results in the collision times tcm(v^) given in Eq.(|^). 
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FIGURE CAPTIONS 



Fig. 1. Analysis of the algebraic decay simulated in pentagon (5-gon) and heptagon 
(7-gon) within the collision flcm space. Symbols - numerical data on the partial weights 
C^'^^'^\(p) and Cj'^'^^\(p) deduced from the observed spectra of the survived orbits N^^'^^\ip) 
and N^j^^\(f) with the help of Eq.(^ and simulated for (f-set orbits with < (f < (fm at 
distinct times t = 20,30rem . Line: the overall-collision-anlle weight Ccm^^ ■ 

Insert left: Points - data Nf^ for the overall survived orbits at late times and their 
regular-orbit analysis with the help of Eq.(|]). 

Insert right: Points - data for (f-set collision time tcslv') simulated within the reduced 
domain < ip < vr/lO. Line - the same predicted in Eq.(|^). 

Fig. 2. Temporal evolution of the survived orbits in rational polygons with small opening 
width (A = 0.05i?) against the reduced time in log-log coordinates. Reduction is given by 
the help of Eq.([I|) for the escape time Xg = 300, chosen common for all cases. Points: 
numerical data for decay of the Qm space phase simulated by A^o = 10^ particles in the 
m-gons (squares) and the correspondent CB (circles). 

Fig. 3. Temporal evolution of the survived orbits in rational polygons with large opening 
width (A = 0.20i?) against the reduced time in log-log coordinates. Notations are the same 
as in Fig. 2. 

Fig. 4. Estimation of the yj-orbit collision time tcm(v^) on the bases of Eq.(A3) for the case 
of m = 3. The regular piece-line orbit a, b, c, d,e... is represented by the infinite straight-line 
trajectory in the triangle LG lattice with the intersection-point sequences 1, 2, 3, 4...n(t, (pkm)- 
The equivalent walls k, the unreduced collision angles fkm, and the axillar angles Qkm are 
shown. 
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TABLE 1 



m 


^(tot) 


^(reg) 


P,(tot) 
J-^ cm 


3 


0.140 


0.116 


0.135 


4 


0.220 


0.219 


0.210 


5 


0.094 


0.086 


0.090 


6 


0.149 


0.139 


0.150 


7 


0.092 


0.069 


0.090 


8 


0.099 


0.096 


0.100 



Table 1. Fitting parameters of the temporal algebraic decay of the collision space ficm 
simulated in the weakly open m-gons with A = 0.05i?. Notations: dm*^ and dm*^^ corre- 
spond to the data on y9-sets observed in the decay spectra C'^^\(p) (see Fig. 1) and averaged 
over, respectively, all collision angles and with excluding singular-orbit angles manifested by 
the high peaks; Dcm*'' = -Dcm^^- the overall-set weights of the algebraic tail given in Eq.(|D 
and derived within the primary relaxation window (see the left insert in Fig.l). 
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